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,  ABSTRACT 

A  combined  Raman-Rayleigh  Udar  (Light  Detection  And  Ranging)  has  recently  been  implemented 
at  the  Air  Force  Geophysics  Laboratory's  ground-based  lidar  station  to  measure  neutral  density  from  the 
lower  stratosphere  to  the  upper  mesosphere. 

Rayleigh  Lidar  reliably  measures  relative  densities  in  the  region  above  30  km.  In  this  region, 
atmospheric  extinction  can  be  neglected  and  backscattering  is  primarily  due  to  Rayleigh  scattering. 
However,  when  density  measurements  are  needed  for  the  lower  stratosphere  two  complications  arise:  the 
contribution  of  aerosol  (Mie  and  Rayleigh)  scattering  to  the  Rayleigh  signal  and  the  effect  of  aerosol 
attenuation.  Vibrational  Raman  scattering,  being  an  elastic  process  for  molecules  only,  can  be  used  to 
resolve  the  first  ambiguity.  The  second  difficulty  requires  an  inversion  technique  to  help  determine  the 
attenuation  profile  from  the  Lidar  signal  and  provide  a  transmission  correction  of  this  signal. 

For  the  lower!  stratosphere,  the  technique  adopted  in  this  laboratory  is  a  three-step  treatment  of 
data  In  step  (1)  Kletr  inversion  is  applied  on  the  elastic  scattering  signal  (Mie  and  Rayleigh)  to  obtain  a 
transmission  altitude  profile.  In  step  (2)  the  molecular  signal  from  the  Raman  Lidar  is  corrected  for 
atmospheric  attenuation.  In  step  (3),  Raman  data  for  below  25  km  Is  spliced  to  Rayleigh  data  lor  above 
25  km  to  give  the  entire  profile  of  neutral  density.  Application  of  this  analysis  to  experimental  data  will  be 
shown  and  discussed. - f  >  f ,  A  . 

\  />  up/ 

1.  INTRODUCTION  ) 


Under  normal  conditions  elastically  backscattered  light  from  the  30-90  km  range  corresponds  to 
that  expected  from  Rayleigh  scattering  by  atmospheric  molecules2.  With  commercially  available  high- 
power  lasers  and  improved  electronic  detection  systems,  Rayleigh  Udar  has  become  a  reliable  tool  in 
determining  neutral  density  in  the  middle  atmosphere.  Furthermore,  with  Increased  laser  power  and 
stronger  signal  strength,  measurements  for  small  intervals  of  signal  integration  are  possible.  This  opens 
up  the  possibility  of  studying  short  time-scale  density  perturbations  such  as  those  caused  by  acoustic 
gravity  waves3. 

Udar  measurements  are  also  difficult  to  calibrate  as  signal  strength  depends  on  transmitting  and 
receiving  efficiencies.  These  vary  with  meteorological  conditions  as  well  as  instrumental  conditions. 
Fortunately,  instrumental  changes  can  be  quantified  and  the  atmosphere  does  not  change  during  the  beam 
round  trip  time.  These  imply  that  after  proper  corrections  for  atmospheric  transmission,  Udar  signals  can 
be  reduced  to  a  relative  neutral  density  profile.  Therefore  only  one  calibration  point  Is  required  for  the 
entire  density  profile.  Slow  changes  of  the  atmosphere  require  independent  measurements  for  calibration. 
Rawinsonde  balloons  can  be  released  so  that  they  floating  above  the  Udar  site  when  attaining  the  altitude 
of  calibration. 

Raman  Udar  has  been  implemented  to  measure  density  profile  for  various  molecular  components 
current  affiliation:  M.  I.  T.  Haystack  Observatory,  Westford,  Massachusetts. 
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of  the  atmosphere4.  Since  Raman  signal  is  spectrally  separate  from  Rayleigh  and  Mie  signals,  it  depends 
only  on  the  density  of  the  molecule  of  interest5.  Together  with  the  Kiett  technique,  this  utility  helps  obtain 
transmission-corrected  absolute  neutral  densities  in  the  regions  below  and  above  the  rawinsonde 
calibration  point.  We  will  discuss  (a)  the  Laboratory  Udar  and  its  operations,  (b)  the  Raman  augmentation 
and  (c)  Rayleigh-Raman  combined  absolute  density  measurements! 

2.  AFGL  LABORATORY  UDAR 
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Figure  1  UDAR  Block  Diagram 

AFGL  has  two  Udar  systems:  a  laboratory  system  (GLEAM,  Ground-based  Udar  Experiments  of 
Atmospheric  Measurements)  and  a  mobile  version  for  field  measurements  (GUNT,  Ground-based  Udar 
INvestigation-Transportable).  The  Raman  capability  has  been  implemented  In  the  GLEAM  system  which 
consists  of  a  Nd:YAG  laser  transmitter  and  a  f/15  36  inch  Dall-Kirkham  telescope.  The  frequency-doubled 
and  -tripled  outputs  of  the  Nd:YAG,  at  S32  and  355  nm  respectively,  is  beam-expanded  with  a  magnification 
of  I5x  and  transmitted  vertically.  The  transmitted  laser  beam  has  a  divergence  of  0.05  mrad  which  is  small 
compared  to  the  telescope  field  of  view.  The  laser  is  operated  at  a  repetition  rate  of  10  Hz  and  pulse 
energies  are  maintained  at  ca.  250  mJ  lor  green  and  60  mJ  for  UV  output  Scattered  radiation  is  collected 
with  a  telescope  field  of  view  which  is  ca.  1.5  mrad.  The  returned  light  is  mechanically  shuttered,  split, 
attenuated  and  coupled  to  a  6  photomultiplier  (PMT)  detection  system.  The  later  consists  ol  a  green  and 
a  UV  channel,  each  channel  having  a  high,  a  medium  and  a  low  altitude  PMT.  The  use  of  beam  splitters 
and  attenuators  ensures  that  the  Lidar  returns  falls  in  the  optimum  dynamic  range  of  each  PMT.  For  most 
applications,  the  return  is  photon  counted.  Photons  are  counted  with  Camac-based  multichannel  scalers 
and  recorded  with  a  personal  computer.  A  safety  radar  is  used  for  aircraft  avoidance.  Fig.  1  shows  the 
Lidar  block  diagram.  For  this  study,  Lie  scaler  dwell  time  is  set  at  2  /re  and  represents  a  range  resolution 
of  0.3  km. 

A 12  inch  telescope  collects  the  Raman  signal.  An  interference  filter  of  10  nm  bandpass  centered 
at  386  nm  is  used  to  spectrally  select  the  first  Stokes  scattered  radiation  from  nitrogen  molecules.  The 
Raman  signal  is  photon-counted  and  integrated  with  a  multichannel  scaler.  Since  the  Raman  signal  is  4 
orders  of  magnitude  smaller  than  the  Rayleigh  signal,  3  orders  for  the  difference  in  scattering  cross 
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sections  and  1  order  for  the  difference  in  telescope  apertures,  only  one  photomultiplier  is  used  for  the 
Raman  return.  The  Rayleigh  signals  from  the  36-inch  telescope  and  Raman  from  the  1 2-inch  telescope 
are  collected  simultaneously.  In  the  basic  data  files,  these  sign  us  are  integrated  over  5  minutes  or  3000 
laser  pulses.  Laser  energies  are  also  recorded  for  post-experiment  power  correction. 

Rawinsonde  balloons  released  at  OOOOZ  in  Albany,  New  York  (130  miles,  West),  Portland,  Maine 
(95  miles,  North  East),  and  Chatham,  Massachusetts  (85  miles,  South  East)  provide  density  measurements 
from  ground  level  to  about  25-30  km.  The  density  profiles  obtained  with  the  slowly  ascending  balloon  are 
path-dependent  and  cannot  be  directly  compared  with  the  real-time  measurements  of  Lidar.  The  calibration 
procedure  ties  the  Lidar  density  to  a  Rawinsonde  data  point  at  25  km. 

3.  TRANSMISSION  CORRECTION  WITH  KLETT  INVERSION: 


To  correct  for  attenuation  in  the  region  below  30  km,  it  is  necessary  to  consider  atmospheric 
attenuation  and  aerosol  (Mie  and  Rayleigh)  scattering.  The  combined  Raman-Rayleigh  Lidar  solves  both 
problems  with  the  following  technique.  We  have  transmitted  radiation  at  two  wavelengths:  532  and  355 
nm.  The  green  signal  is  received  at  the  unshifted  wavelength  for  Rayleigh  data.  The  UV  signal  is  received 
at  unshifted  and  at  Stokes-shifted  wavelengths  -355  and  386  nm  for  nitrogen  molecules,  respectively.  The 
Rayleigh  signal  is  free  of  aerosol  scattering  for  the  region  above  25  km.  Aerosol  anomalies  higher  than 
25  km  such  as  noctilucent  clouds  would  be  detected  in  the  green/UV  ratio  profile.  The  Raman  signal  is 
free  of  aerosol  scattering  but  nevertheless  attenuated  in  the  round-trip  through  the  lower  atmosphere.  The 
transmission  at  355  nm  is  determined  by  the  inversion  of  the  Rayleigh  signal  at  355  nm.  Note  that  the 
round-trip  attenuation  of  the  Raman  signal  is  for  355  nm  wavelength  on  the  way  up  and  386  nm  on  the 
way  down  but  the  difference  due  to  the  wavelength  shift  is  neglected  on  valid  grounds.  Above  25  km,  the 
integrated  transmission  is  expected  to  level  off  and  for  that  reason,  the  532  nm  Rayleigh  signal  is  used  in 
density  analysis  without  correction  for  transmission.  The  transmission-corrected  density  profile  (at  low 
altitudes)  produced  by  Raman  Lidar  is  spliced  to  the  532  nm  Rayleigh  profile  (at  high  altitudes)  to  form  a 
single  profile  from  10  to  80  km.  Figure  2  shows  the  road  map  of  data  analysis. 


The  Klett  technique  is  based  on  the  assumption  that  the  backscattering  coefficient  can  be 
expressed  as  a  function  of  the  extinction  coefficient  In  an  equation,  usually  referred  to  as  the  power  law: 


0  -  const  tfk  (1) 


where  k  is  a  constant  reported  to  be  smaller  than  1  and  greater  than  0.67*.  The  atmospheric  conditions 
of  greatest  interest  to  Lidar  experiments  are  those  where  molecular  scattering  is  stronger  than  aerosol 
scattering  in  the  region  where  density  is  to  be  determined.  Haze  and  optically  thin  cloud  in  the  boundary 
layer  does  not  affect  the  inversion  since  Klett  algorithm  is  worked  downward  and  stopped  at  a  height  of 
10  km.  Under  the  mentioned  conditions,  k  is  expected  to  be  close  to  1  in  the  region  of  interest  Klett 
showed  that  the  extinction  coefficient  along  the  Une-of-sight  depends  on  the  signal  strength  as  described 
by 
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In  the  expression,  S(r)  is  the  range  corrected  logarithmic  photon  count  at  range  r  and  defined  as 


S(r)  -  Infr2 .N(r))  (3) 


where  N(r)  is  the  photon  count  (signal)  at  range  r.  The  use  of  (2)  requires  the  knowledge  of  the  extinction 
coefficient  oM  at  range  rM,  typically  between  25  and  30  km.  in  this  region  molecular  scattering  is 


Figure  2  Data  Block  Diagram 


predominant,  density  is  known  from  Rawinsonde  measurement  and  the  extinction  coefficient  is  readily 
approximated. 

Figure  3  shows  a  typical  355  nm  Udar  unprocessed  data  Laser  energy  is  1 50  mJ /pulse  and  laser 
repetition  rate  is  10  Hz.  Signal  is  received  with  the  36  inch  dia  telescope  and  integrated  for  30  minutea 
The  signal  from  below  8  km  is  modulated  by  the  opening  of  the  mechanical  shutter  and  the  overlap  of  the 
telescope  field-of-view  and  the  beam.  The  shutter  is  fully  open  and  the  overlap  is  total  above  8  km.  Figure 
4  shows  the  corresponding  Klett-inverted  integrated  transmission  profile  at  low  altitudes.  The  values  shown 
correspond  to  round  trip  transmission  between  9.7  km  and  the  point  of  interest .  In  the  inversion,  we  have 
use  0.00027  /km  as  the  extinction  coefficient  at  30  km7.  The  inversion  result  has  been  shown  not  to 
depend  critically  on  this  estimate1.  The  curve  In  Figure  4  indicates  that  density  would  be  overestimated 
by  20%  at  10  km  were  atmospheric  extinction  not  accounted  for  and  the  calibration  tie  made  at  30  km  . 
The  curves  also  indicates  that  the  integrated  transmission  has  levelled  off  by  and  above  25  km.  This  allows 
the  use  of  uncorrected  Rayleigh  signal  above  25  km. 


4.  COMBINED  DENSITY  PROFILES 

Figure  5  shows  the  combined  neutral  density  ratio  to  the  USSA76  standard.  The  profile  is 
measured  at  Hanscom  Air  Force  Base,  on  September  1 9, 1 988  at  0200-0230Z.  The  transmission-corrected 
Raman  density  profile  for  the  10-25  km  range  (lower)  is  spliced  to  the  Rayleigh  10-70  km  density  profile 
(upper).  Both  curves  have  a  1.5  km  vertical  resolution  and  a  30  minute  integration  time.  The  Rayleigh 
profile  for  the  10-25  km  range  is  shown  to  demonstrate  the  effect  of  atmospheric  extinction.  The  valid 
density  profile  is  composed  of  the  Raman  profile  below  25  km  and  the  Rayleigh  profile  above  25  km. 

in  this  laboratory  (result  shown  in  Fig.  5),  Raman  signal  is  obtained  with  the  355  nm  beam  and 
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Figure  5  Combined  Raman-Rayleigh  density  profile.  Raman  profile  (lower  trace  and  extending 
from  10  to  25  km)  has  been  corrected  for  transmission.  Rayleigh  profile  (upper  trace  and  from  10 
to  70  km)  is  also  shown. 
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